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Abstract. When two or more atoms bind to a solid surface, the substrate 
can mediate an interaction between them. In this paper, we use density- 
functional theory to quantify the substrate-mediated pair interaction be- 
tween two adatoms on a compressively strained Ag(lll) surface and on 
unstrained Ag(lll). On the strained surface, the elastic interaction is sig- 
nificant over the short range and leads to a net attraction between two 
adatoms. However, at the longest distances probed, the interaction is pri- 
marily electronic and repulsive. The repulsion can be as strong at 50 meV, 
and it forms a ring-like structure around an atom. On unstrained Ag(lll), 
the interaction is primarily electronic in origin and it is weak relative to the 
interaction found on the strained surface. We calculate the energy barrier 
for an isolated atom to diffuse in each of these systems. For the strained 
surface, the magnitude of the diffusion barrier is comparable to that of 
the adsorbate interaction. We discuss the implications of our findings for 
growth at surfaces. 



1. Introduction 

Adsorbate-adsorbate interactions at surfaces play a key role in controlling 
the structure and properties of interfacial materials. Such interactions can 
be classified as through-space (direct) or through-substrate (indirect), al- 
though a clear distinction is not possible [1]. Quantification of the indirect, 
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substrate-mediated interaction between two adsorbates is the primary sub- 
ject of this paper. 

Substrate-mediated interactions can have an elastic and an electronic 
component. Elastic interactions arise when an adsorbed atom perturbs the 
positions of surrounding surface atoms [2] . Additionally, an adsorbate may 
perturb the charge density near the surface and alter the binding energies of 
neighboring adsorbates - giving rise to an electronic interaction [1, 3, 4, 5, 6]. 
The electronic interaction is oscillatory in nature and, in the limit that the 
separation d between two adsorbates is much greater than the inverse Fermi 
wavevector k^ 1 , the interaction energy AE is predicted [5] to have the form 

AE ~ cos ^ d ) , (1) 
d x ' 

where x = 2 if the substrate has a surface state and x = 5 in the absence 
of a surface state. Oscillatory, charge-density waves associated with the 
electronic interaction have been imaged in scanning-tunneling microscopy 
(STM) studies of noble-metal surfaces [7, 8], which posses a surface state. 
For example, on Ag(lll), kp = 0.08 A [8] and charge density waves can be 
imaged at distances greater than 70 A [8] in STM studies. 

Electronic, substrate-mediated interactions should also be manifested at 
"intermediate" distances between the short separations required for chem- 
ical bonds and the long separations for which Eq. 1 is applicable. For ex- 
ample, Einstein and Schrieffer showed in tight-binding calculations that 
the intermediate-range interaction is present and is, in fact, more signif- 
icant than the long-range interaction [3]. Recent, low-temperature STM 
studies also indicate that this is the case [9] for Cu adatoms on Cu(lll). 
The intermediate range is accessible to modern, electronic-structure calcu- 
lations. In this study, we use density-functional-theory (DFT) calculations 
to quantify the substrate-mediated interaction between two adatoms on a 
(compressively) strained and on an unstrained Ag(lll) surface. We discuss 
the ramifications of these interactions for adsorbate diffusion and growth 
at surfaces. 

2. Model and Methods 

Total energies and forces are calculated using the DFT code described by 
Bockstedte et al. [10], employing fully separable, norm-conserving pseudo- 
potentials of the Kleinman-Bylander form [11]. These were generated using 
the code developed by Fuchs and Scheffler [12]. For our calculations with 
Ag, the s pseudopotential component is taken as local and we use p and 
d projectors. The 4d-electrons are taken as valence states and their wave 
functions are expanded in plane waves, which truncate at a cut-off energy 
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of E cut = 50 Ry for our calculations. We use the generalized gradient ap- 
proximation (GGA) of Perdew, Burke, and Ernzerhof (PBE) [13] for the 
exchange-correlation functional. 

The supercell approach is used to describe the surface. For most results 
presented here, the surface is modeled as a repeating (4 x 2) slab with 16 
atoms per unit cell and a thickness of four layers. The vacuum spacing is five 
interlayer distances, or ~ 12A. Adsorbed atoms are placed on one side of 
the slab and field effects due to the asymmetry are correctly treated in the 
code by a compensating dipole field. For some of the adsorbate interactions 
probed, it is expedient to utilize smaller slabs. Thus, we also use four-layer 
(2 x 2), (4 x 1), and (2 x 1) slabs for some of the calculations. For all of 
the main results presented here, the k-space integration is performed with 
the equivalent of 4 special k points in the full, surface Brillouin zone of 
the (4 x 2) slab. To improve k-space integration, the electronic states are 
occupied according to a Fermi distribution with /cbT=0.1 eV and the total 
energies are extrapolated to zero temperature. 

To verify the accuracy of the pseudo-potential, we calculated the lattice 
constant of bulk Ag in its ground state. For these calculations, we use 
E C ut = 50 Ry and 182 k points in the irreducible portion of the Brillouin 
zone. We find a lattice constant of ao = 4.20A, which is within 1 % of the 
values found by Yu and Scheffler [14] and by Ratsch et al. [15], who both 
used DFT with the GGA functional of Perdew et al. [16]. As an additional 
test of the accuracy, we calculated the interlayer relaxation di2 of the first 
layer of Ag(lll). For these calculations, the first layer atoms are relaxed 
until the forces on them are less than 0.005 eV/A. We find that the first 
layer relaxes inward by 2.0 % relative to the bulk, interlayer spacing db, in 
good agreement with DFT local-density approximation (LDA) results by 
Narasimhan and Scheffler [17], who found an inward relaxation of 1.7%, 
and by Xie, et al. [18], who found an inward relaxation of 1.0 %. Medium- 
energy ion scattering experiments indicate that d\2 contracts by 2.5 % 
relative to db at temperatures below 670 K [19]. A more recent low-energy 
electron diffraction investigation indicates no contraction of d\2 occurs on 
this surface [20]. In previous theoretical studies of Ag(lll) [15, 17, 18], it 
was found that differences between 0(23 and d\, are negligible. Thus, in the 
calculations presented here, ^23 is fixed to d\,. 

We study adsorbate interactions on both Ag(lll) and on a strained 
Ag(lll) surface. For the strained Ag(lll) surface, we employ a bulk lat- 
tice constant of 4.01 A. This value is taken from the results of DFT-GGA 
calculations of the lattice constant of Pt by Ratsch, Seitsonen and Schef- 
fler [15]. We choose this substrate to address experimental issues related 
to Ag epitaxy on a Pt(lll) substrate covered by a monolayer (ML) of 
Ag [15, 21, 22, 23]. To describe the bare, strained Ag(lll) substrate, the 
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top layer of atoms is relaxed and the rest are fixed to the strained, bulk 
positions. This results in an outward relaxation of the top layer by 9.0 
% relative to for the strained substrate. The outward relaxation allows 
for relief of the strain induced by the bulk lattice, which is compressed by 
4.61% relative to bulk Ag, as it is described by our DFT calculations with 
the GGA. 

The total interaction energy AE for a periodic slab containing N adatoms, 
of which M are at binding site a and (TV — M) are at binding site b, is given 
by 

AE = Ef +N - ME a s+1 - (TV - M)E b s+1 + (N - l)E s . (2) 

Here, Eg' +N is the total energy of a slab with N adatoms, E s+1 and E b s+l 
are the total energies of slabs containing one adatom, and Eg is the total 
energy of a bare slab. We obtain Eg from the calculations described above. 
In an attempt to isolate the electronic component of AE, we obtain Eg^ N , 
Eg +1 , and E b s+l as follows: A single atom is placed in a binding site (an fee 
or hep three-fold hollow site) on a relaxed (as described above) bare slab 
and its height is optimized with respect to the fixed substrate. ^From this 
calculation, we obtain Eg +1 and E b s+l . To calculate Eg'^_ N , adatoms are 
placed in binding sites on the relaxed (and fixed) substrate with heights 
fixed to values from the single-adatom calculations. We use this procedure 
to obtain most of the results described below. Additionally, to resolve the 
elastic contribution to AE, we obtain Eg'^ N , Eg +1 , and E h s+l by simulta- 
neously relaxing the adatoms and the first-layer substrate atoms for a few 
trial configurations. 

In the DFT supercell approach, the total interaction energy AE is com- 
prised of interactions between different adatoms in the slab and interactions 
between adatoms in the slab and the periodic-image adatoms. We can ex- 
press AE in terms of these interactions using the lattice-gas Hamiltonian 
approach (see, e.g. [24]), which yields 

AE = 2 X/ V(2) (di,j)<n<rj + v(3 Hdi,j,d iyk , dj^a^k + ... . (3) 

Here, the summations run over all sites i in the slab and all sites j and k 
in the supercell (which includes both the slab and its periodic images), a m 
is unity if site m (m = i,j,k) is occupied and zero, otherwise, d m ^ n is the 
distance between sites m and n, V^ 2 \dij) is the pair interaction between 
two adatoms on sites i and j, and '(dij,di t k,dj t k) is the trio interac- 
tion between three adatoms on sites i, j, and k. We neglect higher-order 
interactions. Another assumption implicit in Eq. 3 is that the interaction 
between adatoms at a fixed distance is independent of whether these atoms 
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Figure 1. Sample lattice-gas Hamiltonian configuration. Atoms in the slab are enclosed 
in the dashed box. The lines with arrows indicated relevant distances between adatoms 
in the slab and periodic-image adatoms. 



occupy fee or hep sites. Actually, in the results presented below, interac- 
tions between two adatoms occupying the same type of binding site are 
always obtained for adatoms on fee sites. We performed one test of this 
assumption for two adatoms on hep sites. This calculation showed that the 
interaction was identical to that found for atoms on fee sites. Finally, in 
calculations designed to resolve the electronic interaction (which comprise 
the bulk of the results presented here) , the adatom binding energies on fee 
and hep sites are virtually equal: The fee site is favored by less than 3 meV. 

To illustrate our approach, we show one of the structures probed in our 
calculations in Fig. 1. In this structure, two adatoms are in the (4 x 2) sur- 
face unit cell at a separation of 3.43 A (distance 2, as can be inferred from 
Fig. 4 below). Each of these adatoms has two, neighboring, periodic- image 
adatoms at a distance of 9.08 A (distance 12 shown on Fig. 4 below). 
From Eq. 3, AE for this configuration is comprised of one pair interaction 
at distance 2 and two at distance 12. Thus, for a given adatom configura- 
tion, we express AE as a sum of pair and trio interactions with unknown 
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coefficients. From 13 different configurations, we obtain a system of linear 
equations and solve these for pair-interaction coefficients up to the 13 th - 
neighbor. We assume that all other interaction coefficients are zero. 

3. Results and Discussion 

Before proceeding to our results, we discuss our tests of their accuracy. 
We test the convergence of our results with respect to E cut , the number 
of k points, the slab thickness, and the exchange-correlation functional for 
different adsorbate configurations on Ag(lll) and strained Ag(lll). The 
adsorbate configurations in these tests are chosen for two major reasons. 
First, for each configuration, the total interaction energy is composed of 
interactions at a single, adsorbate-pair distance. Thus, it is possible to ex- 
tract pair interactions without solution of the full system of equations for 
the lattice-gas Hamiltonian. A second reason for choosing these configura- 
tions is that they posses a single mirror plane - the highest symmetry that 
could be achieved here. The (relatively) high symmetry makes it possible 
to retain up to 16 k points in the full surface Brillouin zone of a (4 x 2) 
unit cell, so that we can assess the convergence of the k-point summations 
for these configurations. For configurations of lower symmetry, 4 k-points 
in the (4 x 2) unit cell is the maximum that we could use with the com- 
putational resources available to us. In Table 1 we show the influence of 
the cut-off energy, the slab thickness, and the number of k points on three, 
pair- interaction coefficients for Ag adatoms on Ag(lll). The distances for 
the pair interactions can be seen below in Fig. 4. For each result in Ta- 
ble 1, each of the three energies necessary to obtain the pair-interaction 
coefficient (cf., Eq. 2) is obtained with the indicated cut-off energy, slab 
thickness, and number of k points. 



TABLE 1. Convergence tests for V (2) (d). 



Distance, d 


iVk 


E cu t, Ry 


Layers 


V^(d), eV 


5 


4 


50 


4 


-0.012 


5 


4 


60 


4 


-0.015 


5 


4 


50 


5 


-0.005 


5 


16 


50 


4 


0.012 


4 


4 


50 


4 


0.004 


4 


16 


50 


4 


0.001 


9 


4 


50 


4 


0.005 


9 


16 


50 


4 


-0.001 
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From Table 1, we see that the tests for distance 5 indicate satisfactory 
convergence with respect to the cut-off energy and the slab thickness. The 
results of our k-point test for distance 5 are less satisfying and there is 
a (relatively) large difference between the pair-interaction coefficients ob- 
tained with 4 and 16 k points. To further probe the convergence of our 
results with respect to the k-point summation, we investigated two ad- 
ditional adsorbate configurations containing distances 4 and 9. For these 
configurations, differences in the pair-interaction coefficients are small and 
indicate a satisfactory k-point convergence. It appears that distance 5 is 
somehow unique. We will discuss this configuration in more detail below. 

To assess the dependence of our results on the exchange-correlation 
functional, we obtain pair-interaction coefficients for three, different adatom 
separations on strained Ag(lll) using the LDA [25]. The pseduopotentials 
for these calculations are generated as described above for the GGA calcula- 
tions, except that the LDA is used for the exchange-correlation functional. 
We test this pseudopotential by calculating the bulk lattice constant for Ag 
using the same conditions as for the GGA calculations. We find ao = 4.06 
A, which is in exact agreement with the value found by Narasimhan and 
Scheffler [17] and within 1% of the value found by Ratsch, et al. [15]. For 
strained Ag(lll), we employ a lattice constant with a value of 3.92 A, as 
found by Ratsch et al. [15] in DFT-LDA calculations of the lattice constant 
for bulk Pt. 

In Fig. 2, we show a comparison of the LDA and GGA results. Here, we 
see a small difference between the pair-interaction coefficients for distances 
4 and 9 and a (relatively) large difference between the LDA and the GGA 
at distance 5. Some of the discrepancy between the LDA and GGA results 
can likely be attributed to differences in strain. The lattice in the GGA 
calculations is compressively strained by 4.61 % relative to the bulk, while 
the strain for the LDA system is 3.21 % (for a reference, the experimental 
strain is 4.2 % [22]). As we will discuss in more detail below, we find that 
pair interactions on Ag(lll) follow the same trend as those on compres- 
sively strained Ag(lll), but they are consistently lower. A comparison of 
the LDA and GGA results on Fig. 2 also indicates this trend. Thus, it is 
likely that strain plays some role in the differences between the LDA and 
the GGA results. 

We now return to our discussion of V^(5), for which we observe large 
discrepancies with different k-point summations (Table 1) and with the 
LDA vs. the GGA (Fig. 2). Although the exact reason for this large error 
bar is unclear, one possible explanation is that trio interactions contribute 
to the total interaction energy of the configuration that we used to ob- 
tain y( 2 )(5). As we see in Fig. 3, this configuration consists of a row of 
atoms occupying alternating fee sites along the (110) direction. A trio in 
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Figure 2. Comparison of LDA and GGA pair-interaction energies V (2) (d) for three 
different adsorbate separations d. 
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Figure 3. Configuration of adatoms used to obtain V^(5). 

this configuration would most likely consist of three adatoms in alternating 
fee sites along the (110) direction. Although we are not able to quantify this 
particular trio interaction, we know that the trio interaction V^ 3 \A, 5,9) is 
attractive by 12 meV [26]. Einstein found in simple, tight-binding calcu- 
lations that the magnitude of the trio interaction is determined primarily 
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by the two, shortest adatom separations in the trio [27]. In light of these 
results, a non-zero trio interaction seems possible for the adatom configu- 
ration in Fig. 3. If this is the case, then one trio interaction is incorporated 
into our pair parameterization for V^ 2 \5) and the combined uncertainties 
in both interactions could lead to the relatively large error bar for this 
interaction coefficient. 

Figure 4 shows the pair interaction for the various distances probed. To 
facilitate resolution of the long-range interactions, the short-range interac- 
tions at distances 1 (the nearest-neighbor fee site to the central adatom) 
and 2 ( the second-closest hep site to the central adatom) are not shown. 
For both surfaces, the interaction is relatively weak at moderate distances 
and it is the strongest and the most repulsive at the longest distances 
probed. For strained Ag(lll) [cf., Fig. 4(a)], the interaction can be as re- 
pulsive as ~ 50 meV at the longest distances and the repulsion forms a 
ring-like structure around the central adatom. The interactions are weaker 
on Ag(lll) [Fig. 4(b)] and do not form a complete ring over the distances 
probed here. Generally, we find that the adatom pair interaction is lower 
for Ag(lll) than it is for strained Ag(lll) - at all distances for which the 
pair interaction is appreciable, the attractions on Ag(lll) are stronger and 
the repulsion is weaker. We note that Bogicevic and co-workers [28] have 
recently obtained interactions similar to ours in DFT-LDA calculations of 
Al on Al(lll) and Cu on Cu(lll). 

Since the surface atoms have the same configuration in all of the results 
discussed above, the elastic interaction due to adsorbate-induced rearrange- 
ment of the surface atoms is suppressed. To assess the role of elastic interac- 
tions, we recompute the pair- interaction coefficients in a few trial structures 
for which both the adsorbate and first-layer surface atoms are fully relaxed 
until the forces are typically less than 0.025 eV/A. 

In Table 2, we show the elastic contribution to the pair-interaction co- 
efficients at several adsorbate separations. This is given as AV^(d), where 
AV^ (d) is the difference between the pair-interaction coefficient obtained 
for a fully relaxed structure and that for the electronic interaction. ^From 
Table 2, we see that the elastic interaction plays a very small role in the 
total interaction energy for Ag(lll). The largest relaxation is for a dimer, a 
weak attraction of 14 meV. The elastic interaction is considerably stronger 
for strained Ag(lll) and is seen to be appreciable up to distance 5 on Fig. 4. 
At distances 4 and 5, the reduction of the pair interaction upon relaxation 
is greater than the electronic interaction [cf., Fig. 4(a)], indicating that the 
elastic interaction is dominant at these distances. At longer distances, the 
elastic interaction is weak and we conclude that the repulsive ring seen on 
strained Ag(lll) is electronic in origin. 

Adsorbate interactions at surfaces can influence both thermodynami- 




Figure 4- Pair interaction as a function of distance from a central adatom (shown in 
black) for strained Ag(lll) (a) and unstrained Ag(lll) (b). 
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TABLE 2. The elastic contribution A V (2) (d) (units of eV) to the pair 
interaction between two adsorbates separated by a distance d. 



Distance, d 


AV (2) (d)- Ag(lll) 


AV (2) (d)- Strained Ag(lll) 


1 


-0.014 


-0.051 


4 


0.003 


-0.016 


5 


0.003 


-0.023 


9 


0.003 


-0.007 


11 


0.003 


-0.003 



cally favored structures and adsorbate diffusion. To assess the impact of 
adsorbate interactions on surface diffusion, we obtain the diffusion barrier 
for an isolated adatom to hop between adjacent three- fold hollow sites on 
both surfaces. In these calculations, both the adatom and the surface atoms 
are relaxed. The diffusion barrier is given by the difference between the to- 
tal energy when the adatom is at a bridge site and when the atom is at a fee 
three-fold hollow site. We find for Ag(lll) that the diffusion barrier is 87 
meV. For strained Ag(lll), we find a barrier of 52 meV. Our barriers are 
in good agreement with those of Ratsch et al. [15], who found 81 meV for 
Ag(lll) and 63 meV for strained Ag(lll) in DFT-LDA calculations. Dif- 
ferences between the barriers that we find and those found by Ratsch can 
be attributed to differences in the exchange-correlation functional (LDA vs. 
the GGA used here). Additionally, there are differences between the lat- 
tice constant predicted for the strained Ag(lll) surface in the two studies. 
As discussed above, the lattice in our GGA calculations is compressively 
strained by 4.61 % relative to the bulk, while the strain for the LDA system 
is 3.21 %. Ratsch et al. found that the value of the diffusion barrier on a 
strained Ag(lll) surface is very sensitive to the degree to which the sur- 
face is strained and that smaller diffusion barriers arise on surfaces under 
greater compression [15]. The difference between our diffusion barrier and 
that found by Ratsch et al. [15] is adherent to this trend. Our diffusion 
barriers are also in good agreement with experimental values [22] for Ag on 
1-ML-Ag/Pt(lll) (60 meV) and on Ag(lll) (97 meV). 

Comparing the diffusion barrier for an adatom on strained Ag(lll) to 
the magnitude of the pair interactions on this surface, we see that the 
two are comparable. Thus, one conclusion arising from this work is that 
substrate-mediated interactions can influence adsorbate diffusion and mor- 
phologies formed during growth at surfaces. Our recent kinetic Monte Carlo 
simulation studies of island nucleation and growth show that this is the 
case [23, 26]. Recent experimental STM studies have begun to probe the 
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role of adsorbate lateral interactions in determining the structure of adsor- 
bate islands at surfaces [9, 29]. Theoretical-experimental interplay in this 
area would be fruitful. 
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